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Abstract 

In this article, we assume that there exists a pseudoscalar DsDso{2317) molecular 
state and study its mass with the molecule-type interpolating current in details using 
the QCD sum rules. The numerical result disfavors identifying the charmonium-like 
CJ ■ state y(4274) as the 75sDso(2317) molecule. 
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^ '■ 1 Introduction 

Ph. 

Ph' In 2009, the CDF collaboration observed a narrow structure (the y(4140)) near the J /ip(t) 

^ , threshold with a statistical significance greater than 3.8(7 in exclusive — ?> J/ij)cj)K~^ 

decays produced in pp collisions [1]. The measured mass and width are (4143.0 ± 2.9 ± 
ly-^ ■ 1.2) MeV and (II.TIj'q ± 3.7) MeV, respectively. There have been several identifications 

J> . for the y(4140), such as the molecular state [2], the tetraquark state [3l H], the hybrid 

^ \ state Oil], the re-scattering effect |7j, etc. The Belle collaboration measured the process 

■ 77 (pJ/i^ for the (j)J/ip invariant mass distributions between the threshold and 5GeV, 
lO ■ and observed no signal for the structure y(4140) (t>J/ip, however, they observed a 
fSJ \ narrow peak (the X(4350)) of 8.8^32 events with an significance of 3.2 o" [8]. The mea- 

■ sured mass and width are (4350. Gt^;? ± 0.7) MeV and {U.^t^^^i ± 4.1) MeV, respectively 
[8]. Recently, the CDF collaboration confirmed the y(4140) in the J/V'</'i^^ de- 
cays with a statistical significance greater than 5(T, the Breit-Wigner mass and width 

, are (4143.4^3:0^0.6) MeV and (l5.3l^°i'^ ± 2.5) MeV, respectively [9], which are con- 

■ sistent with the values from the earlier CDF analysis [1]. Furthermore, the CDF col- 
. laboration observed an evidence for a second structure (the y(4274)) with approximate 

significance of 3.1(7. The measured mass and width are (4274.4l|^ ± 1.9) MeV and 
(32.3l?^J ± 7.6) MeV, respectively [9j. 

In Ref.[10], Liu, Luo and Zhu identify the charmonium-like state y(4274) as the 5- 
wave DsDsQ{2'il7) + h.c. molecule with the spin-parity = 0~, and make prediction 
for the mass of the 5-wave DZ)o(2400) + h.c. molecule as its cousin. In Ref.|llj. He and 
Liu take the y(4274) as the DsDso{22>n) molecular state and study the line shapes of 
the open-charm radiative decays and pionic decays. In Ref.[12], Finazzo, Liu and Nielsen 
study the y(4274) as the L>sl)so(2317) + h.c. molecular state with J^^ = 0"+ using the 
QCD sum rules, and obtain the mass M^^fj^^ = (4.78 it 0.54) GeV. 

The mass is a fundamental parameter in describing a hadron, in order to identify 
the y(4274) as the -Ds-Dso(2317) (or Z)sZ)<jo(2317)) molecular state, we must prove that 
its mass lies in the region (4.2 — 4.3) GeV. The normal threshold of the scattering state 
DsDso{23l7) is about 4287MeV [l3], which is slightly larger than the mass of the y(4274). 
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In this article, we assume that there exists a pseudoscalar DsDsq{2317) molecular state 
in the J invariant mass distribution indeed, and study its mass using the QCD sum 
rules O ITSpl - 

The article is arranged as follows: we derive the QCD sum rules for the mass of the 
y(4274) in section 2; in section 3, numerical results and discussions; section 4 is reserved 
for conclusion. 

2 QCD sum rules for the molecular state y(4274) 

In the following, we write down the two-point correlation function H{p) in the QCD sum 
rules, 

n(p) = iy"d^2;e^f^(0|r{j(x)jt(0)} |0), (1) 
J{x) = c{x)i^^s{x)s{x)c{x) , (2) 

we choose the pseudoscalar current J{x) to interpolate the molecular state y(4274). There 
are two additional pseudoscalar currents ri^{x) and J5(x), 

r/5(x) = c(a;)7^75s(x)s(x)7^c(x) , 

J5{x) = c{x)(j^^-f5s{x)s{x)a''''c{x) , (3) 

which interpolate the pseudoscalar i:>*(2112)Z)*i(2460) and L>^i(2536)l)si(2536) molec- 
ular states, respectively, here we have neglected the mixing between the Z)*]^(2460) and 
Z)si(2536) [16], and take them as the ^Pi and ^Pi states respectively in the non-relativistic 
quark model considering the masses M£,* k, M£)^^. The normal thresholds of the scat- 
tering states DsDsQ{2m), L>*(2112)Z)*/(2460) and L>si(2536)l),i (2536) are 4287 MeV, 
4572 MeV and 5070 MeV, respectively [13]. It is impossible for the L>,i(2536)L>,i(2536) 
molecular state has the binding energy —796 MeV so as to reproduce the mass of the 
F(4274). In Ref.[6], we perform detailed studies of the L»*(2112)i5*(2112) molecular state 
using the interpolating current, 

r/(x) = c(j;)7^s(x)s(j;)7^c(x) , (4) 

and obtain the value M^*^* = (4.48 ± 0.17) GeV or (4.43 ± 0.16) GeV, which is larger 
than the normal threshold of the scattering state Z)*(2112)l)*(2112), and draw the con- 
clusion tentatively that the (2112)1)* (2112) is probably virtual state not related with 
the y(4140). Compared with the scalar current r/(x), the pseudoscalar current r]^{x) has 
an additional Dirac matrix 75, which leads to the terms m^{ss)'^, m'^{sgsCrGs)'^ , and some 
terms mdss), mdsgsCrGs) change their sign, and the convergent behavior in the operator 
product expansion becomes worse, we have to choose much larger Borel parameter and 
postpone the threshold parameter to much larger value, and obtain the value M^),^). , 
which is much larger than the M^,*^*. It is difficult to reproduce the mass of the y(4274) 
as the pseudoscalar D* (2112)1)*]^ (2460) molecular state using the QCD sum rules. 

^In preparing the article, Ref.[12) appears. We take into account the vacuum condensates adding up 
to dimension-10 consistently, while in the second version of Ref. [12], some vacuum condensates such as 
(ss){^i^), (sgsaGs}{^^^) are still neglected. 
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We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operator J{x) into the correlation function Il{p) to obtain the 
hadronic representation \14:\ [T5] . After isolating the ground state contribution from the 
pole term, which is supposed to be the y(4274), we get the following result, 

where the pole residue (or coupling) Ay is defined by 

Xy = (o|j(o)|y(p)). (6) 

The contributions from the two-particle and many-particle reducible states are small 
enough to be neglected [IT]. The scattering state l)sDso(2317) has the same quantum 
numbers as the current operator J(x), the corresponding two-particle reducible contribu- 
tion can be written as 



n2(p) = ^%^D.ol 7^A- ..r/ ^ + (7) 



(27r) |^g2 _ ^2^J (j) _ qY _ 
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where the pole residue (or coupling) Xq^jj^^ is defined by {0\J{0)\DsDsq{p)) = Xjj^d^^^- In 
the soft Dg meson limit, the pole residue can be estimated as 

{0\jmDsDso{p)) = -^{0\[Q5,Jm\Dso{p)), 

IDs 

Q5 = j (i^2;s(x)7o75c(x) , (8) 

the is the decay constant of the Dg meson. We can carry out the commutator in 
Eq.(8) and obtain a current consists of four quarks. The coupling of the scalar meson 
-Dso(2370) with a four-quark current should be very small [n\ llSj. In this article, we take 
the -Dso(2370) as the conventional two-quark meson. 

The J{x) maybe also have nonvanishing couplings with the following two-particle scat- 
tering states, such as the Xc0?/(548), Xc0^'(958), J/ip7]{5A8), J/V'r/'(958), J/V'/( (1510), 
/ic/ii(1380), Xci</'(1020), r]c4>{1020), 77c/q(1500). The corresponding normal thresholds 
are M^^^r, = 3962.6 MeV, M^^^^/ = 4372.53 MeV, Mj/^^ = 3644.77 MeV, Mj/^^/ = 
4054.7 MeV, Afj/^j/ = 4614.92 MeV, Mh^hi = 4911.42 MeV, M^^^^ = 4530.12 MeV, 
M^^0 = 3999.76 MeV and M^^j/ = 4485.3 MeV, respectively [H]. We can write the 
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interpolating current J{x) in the following form with the Fierz re-ordering, 

1 _ _ i _ 

J{x) = -— c(x)c(2:)s(x)z75s(a;) - — c(a;)7ac(a;)s(a;)7"75s(a;) 

-^c(a;)(7a/3c(a;)s(x)(7"^75s(x) + ^c{x)-falbc{x)s{x)'y°' s{x) 

1 1 A" A" 

-z^c{x)i'yr,c{x)s{x)s{x) - -c{x)—c{x)s{x)—i^^s{x) 

i A'' A'* i A" A" 

-lf{x)^lac{x)s{x)—-i'^-i^s{x) - -c(x)ycj„^c(x)s(x)ya"^75s(x) 

i X"" A" 1 A'' A" 

+if{^-)^lalbc{,x)s{x)—^°'s{x) - -c{x)—i-i5c{x)s{x)—s{x), 

= -^4p{x) - ^4a{^) - y^4t{x) + Y^J%{x) - ^^J%{x) 

-\4p{x) - \Jva{x) - \Jtt{x) + \Jav{x) - \jps{x) , (9) 

where the and are the color-singlet and color-octet currents, respectively. Here we 
have used the following identity, 

1 / A'* \ / X"'\ 

to perform the re-arrangement in the color space. There is another identity, 

m.«-Kt).(t) ,-I(t),,(t) ■ 

to express the current J{x) into a series of color-octet currents, the identities in Eqs.(10- 
11) can be changed into each other. The couplings to the lowest scattering states below 
the normal i)sL'so(2317) threshold can be estimated as 

{0\JlpmXcor]{p)) = -f(Ol[Q5,J^p(0)]|Xco(p)), 

(0|4^(0)|J/^r?/r?'(p)) = -^(0|[Q5,4a(0)]|J/V'(p)), 

Jv/v' 

{o\j%vmvcm) = -^(oi[q5,jV(o)]|0(p)), (12) 

in the soft pseudoscalar mesons limit, where the Q5 are the axial-charges in the due 
channels, the f^f, frj' and /^^^ arc the decay constants of the pseudoscalar mesons. If we 
carry out the calculations for the commutators in Eq.(12), we obtain the currents consist 
of four quarks. The couplings of the four-quark currents to the two-quark mesons are 
supposed to be small, the contaminations from the two-particle reducible scattering states 
can be neglected. 

On the other hand, the pseudoscalar charmonia rjc, rf^, rj", ■ ■ ■ have Fock states with 
additional qq components besides the cc components. The current J(x) maybe have non- 
vanishing couplings with the pseudoscalar charmonia, those couplings are supposed to be 
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small, as the dominating Fock states of the pseudoscalar charmonia are the cc components. 
The contaminations from the conventional charmonia can also be neglected. 

The molecular current J{x) can also be expressed in terms of a series of diquark- 
antidiquark type currents with the rearrangement in the Dirac spinor space and color 
space, in other words, the c{x)Tc{x)s{x)T' s{x) type currents in the J{x) can be recasted 
into a special superposition of the diquark-antidiquark type currents themselves. The 
molecular current c{x)Tc{x)s{x)T's{x) maybe have non-vanishing couplings with a scries 
of tetraquark states, we cannot distinguish those contributions to study them exclusively, 
and assume that the current c{x)Tc{x)s{x)r' s{x) couples to a particular resonance, the 
molecular state (for example, the XcO^^ J/i^V^ Vc^Py etc), which is a special superposition of 
the tetraquark states (or has the tetraquark states as its Fock components), and embodies 
the net effects. On the other hand, one can resort to the tetraquark scenario instead of the 
molecule scenario to study the possible contaminations. However, it is a very hard work, 
as the predicted masses of the tetraquark states from different theoretical approaches differ 
from each other greatly, and the tetraquark states have not been confirmed experimentally. 

In the following, we briefly outline the operator product expansion for the correlation 
function Il{p) in perturbative QCD. We contract the quark fields in the correlation function 
n(p) with Wick theorem, obtain the result: 

U{p) = i J d^xe'P '^Tr [i^5Sab{x)i-i^CU-x)] Tr [Cmn{x)Snm{-x)] , (13) 

where the a, b, m and n are color indexes, and substitute the full s and c quark propagators 
Sab{x) and Cab{x) into the correlation function Ii{p) and complete the integral in the 
coordinate space, then integrate over the variables in the momentum space, and obtain 
the correlation function IL{p) at the level of the quark-gluon degrees of freedom. 

Once analytical results are obtained, then we can take the quark-hadron duality and 
perform Borel transform with respect to the variable = —p^, finally we obtain the 
following sum rule: 

_a4 rso _ , 

Aye aP" = / dsp{s)e , (14) 

where 

A / \ I B ' lOLfiGG 



P{s) = Po{s) + p{ss) (s) + P(ss)2 (s) + P(GG> («) + P(GG) («) ( 



TT 



+P(GGG){s){gtGGG) , (15) 

the lengthy expressions of the spectral densities po{s), P{ss}{-'^): P{ss)'^{^)^ PfGG)^^^' P^GG)^^^ 
and P(^GGG)i'^) presented in the appendix. In this article, we carry out the operator 
product expansion to the vacuum condensates adding up to dimension-10 and take the 
assumption of vacuum saturation for the high dimension vacuum condensates. 

Differentiate Eq.(14) with respect to -jg^-, then eliminate the pole residue Ay, we can 
obtain a sum rule for the mass of the F(4274), 

f/i? , y> ds ,'L .r2\ Pis)e~M^ 
^2 ^ J4(r». + ,»,,)- di-i/M^)Py^ >^ _ ^^g^ 
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3 Numerical results and discussions 



The input parameters are taken to be the standard values (qq) = —(0.24 ± 0.01 GeV)^, 
{ss) = (0.8±0.2)(gg), {sgsaGs) = mg(ss), = (0.8 ± 0.2) GeV^, (^^) = (0.33 GeV^ 
{g^GGG) = 0.045 GeV^ = (0.14±0.01) GeV and nic = (1.35±0.10) GeV at the energy 
scale /u = 1 GeV O EH [HI [20] . 

In the conventional QCD sum rules [14^ I15j . there are two criteria (pole dominance 
and convergence of the operator product expansion) for choosing the Borel parameter 
and threshold parameter sq. We impose the two criteria on the DsDsq{2317) molecular 
state to choose the Borel parameter and threshold parameter sq. 

In Fig.l, we plot the contributions from different terms in the operator product expan- 
sion. From the figure, we can see that the contributions change quickly with variations of 
the Borel parameter at the region < 3 GeV^ , which does not warrant a platform for 
the mass, see Fig. 2. In this article, we can take the value > 3 GeV^ tentatively, and 
the convergent behavior in the operator product expansion is very good. 

In Fig. 3, we plot the contribution from the pole term with variations of the threshold 
parameter sq. From the figure, we can see that the value sq < 27GeV^ is too small to 
satisfy the pole dominance condition. If we take the values sq = (28 — 30) GeV^ and 
= (3.0 — 3.6) GeV^, the pole contribution is about (44 — 68)%, the pole dominance 
condition is well satisfied. The Borel window changes with variations of the threshold 
parameter sq, in this article, the Borel window is taken as 0.6 GeV^, which is small enough. 
If we take larger threshold parameter, the Borel window is larger and the resulting mass 
is larger, see Fig. 2. In this article, we intend to obtain the possibly lowest mass which 
is supposed to be the ground state mass by imposing the two criteria of the QCD sum 
rules. In the Borel window = (3.0 - 3.6) GeV^, the main contributions come from 
the perturbative term + the {ss) term + the {sgsCrGs) term, while the dominating one is 
the perturbative term, the convergent behavior in the operator product expansion is very 
good, see Fig.l. 

Taking into account all uncertainties of the input parameters, finally we obtain the 
values of the mass and pole residue of the y(4274), which are shown in Figs. 4-5, 

My = 5.08l[!;^^ GeV , 

Ay = 5.68tlil X 10"^ GeV^ . (17) 

From Fig. 5, we can see that at the value = 3GeV^, the uncertainty of the pole 
residue is too large, the platform is not flat enough, we can take a smaller Borel window, 
M2 = (3.2 -3.6)GeV2, then 

My = BmtoilGeY, 

Ay = 5.56l?:°2 X 10-^ GeV^ . (18) 

In the QCD sum rules, the high dimension vacuum condensates are always factorized 
to lower condensates with vacuum saturation, factorization works well in large Nc limit. In 
the real world, Nc = 3, there are deviations from the factorable formula, we can introduce 
a factor k to parameterize the deviations, 

(ss)^ , {ss){sgsCrGs) , (sgsCrGs)'^ — > k(ss)^ , K{ss){sgsaGs) , K{sgsCrGs)^ . (19) 
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Figure 1: The contributions from the different terms with variations of the Borel pa- 
rameter in the operator product expansion. The A, B, C, D, E and F corre- 
spond to the contributions from the perturbative term, {ss) + {sQsCtGs) term, ( "°^'^ ) 
+ (^^) [{ss) + {sgsoGs) + (ss)^] + {glGGG) term, (ss)^ term, {ss){sgs(TGs) term and 
(sgsCrGs)'^ term, respectively. The notations a, /3, 7, A, r and ^ correspond to the threshold 
parameters sq = 26 GeV^ 27 GeV^, 28 GeV^, 29 GeV^ 30 GeV^ and 31 GeV^, respectively. 
Here we take the central values of the input parameters. 
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Figure 2: The mass with variations of the Borel parameter M"^ and threshold parameter 
So- The notations a, /3, 7, A, r and ^ correspond to the threshold parameters sq = 26 GeV^, 
27GeV2, 28GeV^ 29GeV^ SOGeV^ and 31GeV^ respectively. Here we take the central 
values of the input parameters. 
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Figure 3: The contribution from the pole term with variations of the Borel parameter 
and threshold parameter sq- The notations a, /3, 7, A, r and ^ correspond to the 
threshold parameters sq = 26GeV^ 27GeV^ 28GeV^ 29GeV^ SOGeV^ and SlGeV^, 
respectively. Here we take the central values of the input parameters. 
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In Fig. 6, we plot the mass My with variations of the parameter k at the interval 
K = — 2. Prom the figure, we can see that the value of the My changes quickly at 
the region M2 < 3.0 GeV^, and increases with the k monotonously. At the interval 
= (3.2 - 3.6) GeV^, the value k = 1±1 leads to an uncertainty about lo;^7GeV, 
which is too small to smear the discrepancy between the present prediction (5.06 GeV) 
and the experimental data (4.274 GeV) |9j. In the QCD sum rules for the masses of the 
p meson and nucleon, k > 1 [21J. If the same value holds for the molecular states, the 
deviation from the factorable formula means even larger discrepancy between the present 
prediction and the experimental data. 

The central value My = 5.06 GeV is about 770 MeV above the DsDso{23l7) threshold 
[13j . the DsDso(2317) is probably a virtual state and not related to the charmonium-like 
state y(4274). The present prediction is considerably larger than the value M^^^,^^ = 
(4.78 ± 0.54) GeV p]. In Ref.[l2], Finazzo, Liu and Nielsen take the threshold param- 
eter as So = {My + 0.5 GeV)^, and adjust the threshold parameter and Borel parameter 
to obtain the Borel window and reproduce the relation sq = {My + (0.4 ~ 0.6) GeV)^. 
While in the present work, we search for the threshold parameter and Borel parameter by 
imposing the two criteria of the QCD sum rules, and try to obtain the ground state mass 
(which is not necessarily the same as My) in the Borel window, as the tetraquark states, 
irrespective of the molecule type and the diquark-antidiquark type, have not been firmly 
established yet, and we have no confidence to take the ground state as the y(4274). 

The flux-tube model [22], the lattice QCD [23j and the QCD string model [24J predict 
that the masses of the low lying hybrid charmonia are about (4.0 — 4.2) GeV, (4.0 — 
4.4) GeV and (4.2 — 4.5) GeV, respectively, which are consistent with the experimental 
data. While the QCD sum rules indicate that the masses of the ground-state hybrid 
charmonia with J^'-^ = O"*""*", and 1"' are 5.4 GeV, 5.8 GeV and 4.3 GeV, respectively 
[25j . which disfavors identifying the y(4274) as the O"*"*" hybrid charmonium because the 
CDF collaboration fitted the experimental data to an 5*- wave Breit-Wigner resonance [9]. 

The hybrid mesons usually decay to an S'-wave and a P-wave meson pair, and the 
couplings to two S-wave mesons are suppressed [26]. If the y(4274) is a hybrid char- 
monium, the decay y(4274) — >■ J can take place through the final-state re-scattering 
mechanism, y(4274) DsDso{2'ill) — > J/ipcp, and the decay to two photons should be 
forbidden or very small [27], which is in contrary to the X{4350) observed by the Belle 
collaboration in the (j)J/ "0 invariant mass distributions in the process 77 — >■ (pJ/ ip ^ . The 
X(4350) has been tentatively identified as the 2+"'" cscs tetraquark state [31 [8], the (or not 
the) D*Dso molecular state ([25]) [51[2U1[3U] , the P-wave charmonium Xc2 [31]' scalar 
cscs tetraquark state (321 [33]) the scalar cc — D^D* mixing state [32], etc. 

In Ref.[34], we study the mass spectrum of the scalar hidden charmed and hidden 
bottom tetraquark states in a systematic way using the QCD sum rules, and observe that 
the scalar-scalar type and axial-vector-axial-vector type scalar cscs tetraquark states have 
the lowest masses, about (4.44 it 0.16) GeV, while the pseudoscalar cscs tetraquark states 
have much larger masses. The y(4274) may also be an scalar cscs tetraquark state, and 
decay y(4274) — > J/^<j3 takes place through fall-apart mechanism with rearrangement in 
the color space. More experimental data are still needed to identify the new charmonium- 
like states. 
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Figure 4: The mass with variations of the Borel parameter 
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Figure 5: The pole residue with variations of the Borel parameter 
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Figure 6: The mass with variations of the parameters k and M^, other parameters are 
taken to be the central values. 
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4 Conclusion 



In this article, we assume that there exists a pseudoscalar 5<jD<jo(2317) (or DgDsoi'^'i^^)) 
molecular state in the J/ip(p invariant mass distribution indeed, and study its mass with 
the molecule-type interpolating current in details using the QCD sum rules. The numerical 
result indicates that the mass is about My = 5.061q"22 GeV, which is inconsistent with the 
experimental data. The DsDso{2317) (or DsDso{2317)) is probably a virtual state and not 
related to the charmonium-like state y(4274), and we cannot draw a definite conclusion 
with the QCD sum rules alone. 
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Appendix 

The spectral densities at the level of the quark-gluon degrees of freedom: 
3 



Po{s) = 



20487r6 



raf rl-a 

da I 

Jai J Pi 



dPaPil -a-py-'is- mlf{ls^ - Qsm^ + m^) , (20) 



P(ss){s) 



3ms{ss) f"f 
647r4 



paf pi— a 

/ da dpaf3{l - a- P){10s'^ -12sml + 3fhi) 

Jai J Pi 



1287r4 



327r4 



da / d/3a/3 [6(2s - m^) + s'^S{s - m^)] 

ai J Pi 

■a; ,-l-a ^ 

da I dp{s — m^) 



P^ 



1287r4 



/ da 

Jai J Pi 



,fl~2 3msml{sgsaGs) 



Jai 



da, (21) 



m1{ss)^ 
167r2 



J a 



647r2M2 



Jai 



da 



s — 



f 

Jai 



da 



2M2 
1 



+ 



a{l — a)s^ 
4M4 



a{l-a) M2_ 



+ -n^ o{s — m^) 



(22) 
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5127r4 



n da r f ^ + ^ ) (1 - a - 

Jai J Bi 



/32 a2 



2s — + -^-i^ls — m; 



+ 



10247r4 



da 



l-a 



dp{a + -a- pfilOs"^ - 12sml + 3m^) 



msm1{ss) 
1927r2 



da/ d/3f ^ + ^ ) (l-a-/3) 



msml{sgsaG s) ["f 
23047r2M6 

ms'm'^{ss) 
1927r2M2 



ft V 

da f d(3 ' ^ 

hi 



c2 n 



1 + 



M2 2M4 



5(s - m^) 



da 



l-a 



dl3 



ft 



msm1{ss) 



647r2 

ms{ss) 



da 



fi J /" 1 — a 

da / d/3(a + /3) 

a, J ft 

da 



1287r2 



3 + 2s + 



2M2 



7687r2M4 
uism^isgsaGs) f"f 



2567r2M2 



d/3 

da / d/S 

a-i J Pi 



1 1 

+ 



a/33 ^3^ 



J(s-m^) 



s = 



288M4 



da 



1 



a/32 Q,2^ / \ 
1 



a^ (1 — a)3 



m; 



7687r2M4 
^(ss)2 r"/ 



1 



96M2 



da 



1 1 

a2 (1 — a 



msml{sgs(jGs) 
2567r2M2 

ms{sgsaGs) ["'f 



da 



1 

a^ (1 — a)3 

(5(s - mj) 
1 



5(s -mj 



1 



+ 



7687r2 



da 



a2 (1 — a)2 
2 n 



(5(s - m^) 



s s 
1 + ^ + 



M2 2M4 



(23) 



(24) 



^2 raf rl-a 



d/3 

■l-a 



a /^ 

+ ^ 



{l-a-pf 



1+ ^ + 



3 6M2 



(5(s - m^) 



1 I ,o('^ 



81927r6 

2 

"40967r6 



da 



l-a 



ra^ nl—a 

I da I 

J ai J ft 



d/3{^ + ^)(l-a-/3)^ 



2s — + -^(^(s — 



6s — 3mg + — 5(s — m; 



3 + 2s + 



2M2 



(5(s-mJ , 
(25) 
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wnere a/- = — ^^-^ — ai = — ^—^ — Pi = rrir = - — ^ — -, m„ = n \ - 
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